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TRANSIENT PRESSURE RESPONSE IN A FLUID SYSTEM 
AS A VIBRATING VALVE' CLOSES 


By 

Cecil A. Ponder, Jr. 

George C. Marshall Space Flight Center 
Huntsville, Alabama 
ABSTRACT 

Equations are derived that describe the unsteady flow in a pipe- 
line when a vibrating valve closes; a pressure oscillation occurs in 
addition to .water hammer. The results of analytical and experimental 
investigations show that the pressure oscillation causes an increase 
in the maximum transient pre s sure- during valve closure that is pro- 
portional to the frequency and amplitude of valve vibration. • 

Because the boundary conditions and natural frequencies may vary, 
each fluid system must be investigated to determine 'the relationship's 
between the pressure oscillation and the frequency and amplitude of 
the valve vibration. The analytical methods reported can be used for 
such an investigation. 
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TRANSIENT* PRESSURE RESPONSE IN A FLUID SYSTEM 
AS A VIBRATING VALVE CLOSES 

By Cecil A. Ponder, Jr. * 

George C. Marshall Space Flight Center 

SUMMARY 

The transient pressure response in a pipeline caused by the 
closure of a vibrating valve was investigated in order to study the 
effects of the vibration. The vibrating valve causes a pressure oscil- 
lation which increases the maximum transient pressure. Analytical 
expressions describing the transient flow were derived, and the bound- 
ary conditions of a test facility were used in a solution of the equations . 

In an experimental investigation, transient pressures were recorded 
as the valve closed while vibrating at several different frequencies and 
amplitudes, including the natural frequency of the system. The time of 
valve closure was varied between . 250 and 1. 00 second. Two suction line 
configurations were used. One was a uniform suction line of 4-inch 
diameter; the other had a 6 -inch diameter to 4-inch diameter concentric 
reducer midway between the reservoir and valve. 

The results of the investigation show that the maximum transient 
pressure is increased only slightly by low level -vibration, but as the 
vibration increases, the-re is a significant increase in the maximum 
transient pressure. The acceleration at which this increase in the 
maximum transient pressure becomes significant depends upon the 
natural frequency of the system, and the .frequency and amplitude of the 
forced vibrations. With the valve completely closed, a large pressure 
oscillation was obtained during vibration. 

INTRODUCTION 

Aerospace engineers have recently begun to study the effects of 
fluid vibrations on space vehicles. Because high pressures require 
more rigid ducts and structures, there is particular interest in the 
magnitude of pressure oscillations produced by the fluid vibrations. 

The transient pressure that is caused by the closing of a valve has 
historically been called water hammer. When the closing valve is also 
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vibrating, the water hammer couples with pressure oscillations, and 
a higher maximum pressure results. These pressure oscillations, 
if transmitted throughout the fluid system, create structural problems. 
Because there was a need for understanding the coupling of fluid- 
vibrations and water hammer, the transient pressure response of 
the closure of a vibrating valve was investigated (Ref . 1 and 2). 
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HISTORY OF WATER HAMMER INVESTIGATIONS 


There have been numerous analyses of water hammer. A discus- 
sion of various methods of reducing water hammer by Michaud (1878) 
was the first investigation of the phenomena. In 1897, N. E. Joukowsky 
first connected water hammer with acoustic wave action (Ref. 3 and 4). 
His analysis was based on instantaneous valve closure and a uniform 
pipeline through which an inviscid fluid was* flowing. He derived the 
equation: 



AV 


(1) 


Joukowsky conducted an extensive program to verify his equation. 
However, the effect of compound and branched pipes, or valve closures 
of finite time was not included in his analysis. 

At about the same time Lorenzo Allievi, an Italian engineer, made 
a thorough study of the water hammer caused by linear valve closures 
of finite time (Ref. 5). He constructed a chart from which the maxi- 
mum surge pressure during a linear valve closure could be obtained. 

One of the first water hammer studies in the United States was 
.conducted by N. R. Gibson (Ref. 6). He used the method of arithmetic 
integration to give the same results that Allievi- had obtained, yet he 
had no knowledge of Allievi' s studies. E. E. Halmos pointed out the 
identical results, and in 1925, he published an English translation of 
Alleivi's works. 

Several other investigations of water* hammer were made-, and in 
1927, Ray S. Quick published a comparison of the various theories 
(Ref. 7). These theories, however, were merely amplifications of 
the basic ideas of Joukowsky and Allievi. 

Nonuniform gate closure was investigated in 1928, by S. Logan 
Kerr (Ref. 8) in connection with his investigations of valve closures 
from partial openings. R. W, Angus extended Joukowsky' s theory to 
include a thorough study of water hammer in compound and branched 
pipe systems (Ref. 9). He relied heavily on graphical solutions. 

J. Parmakian published Waterhammer Analysis (Ref. 10) in 1955. 
It included a derivation of the basic linearized water-hammer equations 
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with an abrupt pressure loss at the end of the pipeline to account for 
the frictional losses. Graphical methods used to solve water-hammer 
problems were thoroughly investigated by Louis Bergeron (Ref. 11). 

F. M. Wood used Heavisides' operational calculus as an analytical 
method of predicting the transient pressures in the pipeline (Ref. 12). 
Although the nonlinear frictional loss term of the water-hammer 
equations was linearized, he neglected the other nonlinear terms of 
the wave- equations. Wood's analysis gave the results in terms of only 
the pressure and velocity surges that were produced. 

R. George Rich improved Wood's analysis by using the La Place- 
Mellin transformation so that he worked directly with the total pressure 
in the pipeline (Ref. 13). He also used the linearized friction loss and 
ignored the other nonlinear terms of the water-hammer equations. 

Many authors, including Parmakian, Wood and Rich believed that the 
solution of the water-hammer equations with the nonlinear terms was 
impossible (Ref. 14), 

However, in I960, M. Lister published an article on the solution 
of hyperbolic partial differential equations (Ref. 15). The article 
showed that the theory of characteristics could be used to solve partial 
differential equations with two dependent and two independent variables. 

V. L. Streeter and C. Lai used this theory of characteristics with 
an approximation to permit calculation of the transient pressure at 
specific time intervals (Ref. 16). 

The solution used in this investigation is based upon the theory of 
characteristics as suggested by M, Lister. A grid of characteristics 
was used, however, to eliminate errors introduced by the method of 
specific time intervals. The spacing of time was controlled by the 
number of characteristic grid lines used. 

By using the theory of characteristics to solve the wave equations, 
the transient pressure and fluid velocity are calculated not only at the 
ends of the pipeline, but also along the length of the pipe. There is 
also no need for linearizing or neglecting the nonlinear terms. Because 
the primary purpose of this investigation was to study the effect of a 
vibrating valve, this additional boundary condition was included in the 
analysis. 
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ANALYSIS 


Unsteady Flow in a Uniform Pipeline 


The solution of the transient pressure is based on a one-dimen- 
sional analysis of the transient flow through a uniform pipeline, with 
boundary conditions given at each end of the pipe, as well as the initial 
steady state conditions along it. The fluid system investigated is shown 
in FIG 1. It has a uniform pipeline with a reservoir at one- end and a 
valve at the other end. The basic equations describing the flow through 
the pipeline are: 

a. The momentum equation, Ref. 14,. (see Appendix A for ’deriv- 
ation) 

_av_ V v_»v + J2L W+ jP = 0 (2) 

8t 3x 2D P 9x 

b. The equation of continuity for unsteady flow of a compressible 
fluid (Ref. 10) 


dp 

at 


+ p 


av 

ax 


+ v 


|£. = 

8x 


(3} 


c. The equation of state for a compressible liquid (Ref. 16) 


8P J3_ (4) 

ap p 


W-here j3 is the effective bulk modulus of the liquid in an elastic 
container defined as (Ref; 10). 


where values 
C, = 5/4 - [i 


P= i 

_L + c i p ' 

K eE 

of C p are as follows: 

for a pipe anchored at the upper end without 
expansion joints 


(5) 


5 



C ( = 1 - jj. 2 . for a pipe anchored against longitudinal movement 
throughout its length 

C 1 = 1 - p/2 for a pipe with expansion joints 

It should be noted that P is the piezometric pressure, and there- 
fore, includes pressure variations caused by elevation; that is, 

P = pz + p, where p is the local static pressure and z is the elevation 
above a datum plane. 

Combining the equation of continuity (3) and the equation of state 
(4) 


9P 

8t 


+ 


pH. 

ax 



( 6 ) 


This equation, along with the momentum equation (2), forms the 
basic wave equations used in the analytical study of water hammer. 

By using the theory of characteristics, a set of characteristic 
equations is derived in Appendix B. These characteristic equations 
define two sets of curves as shown in FIG 2. One set of curves has 
a positive slope and the other set has a negative slope (identified as 
C + and C"). 

The characteristic equations are: 

Along C + 


( x o - X A> = < V A + a A> <*0 - f A> ( 7 ) 


(V ° ' VA) + "fe <P ° ' ?A) - *A)-0 




(8) 

Along C 

< x o - X B> = ( V B - a B> (to - t B ) 

(9) 

(▼o - 

Vb) " P B a B (P °- P B) + — B 2 j VB| (to-t B )=0 

(10) 
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With the fluid velocity (V), pressure (P), location (x), time (t), 
and density (. p) given at points A and B, equations (7), (8), (9), and 
(10) can be solved simultaneously for the four unknown quantities 
(P, V, x, t)' at point o. It is also necessary to determine the density 
at point o. The equation of state (4) can be integrated from point A 
to point o. to give the density. 



Po = ^A e 

by using the series expansion (Ref. 17) 

5c 2 tt 3 Y n. 

e x =l+x+2r + 3!~ + •■■ + ~ U 4 > 



the density at point o can be approximated with a high degree of 
accuracy as 



It should be noted that ^ — Q ■ ~ ?A..j is of the order of 10"^ so that 
higher order terms may be neglected. As shown in the FIG 2, the 
quantities at point o are later used to determine the unknown quantities 
at point S. The number of points for which solutions can be obtained is 
one less than the number of points for which the conditions initially are 
known. This problem can be eliminated, however, if the boundary 
conditions are known. * 


Boundary Conditions at the Reservoir 


It is assumed that the reservoir used is of sufficient capacity to 
provide a constant total pressure at the intersection of the pipeline 
and the reservoir. Therefore, returning to FIG 2 

■ V T r 


-•total 


= P L+ 2g~ 


"i/L 


CONSTANT 


( 16 ) 
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This boundary condition is used in conjunction with the equation for 
the C~ characteristic to evaluate the unknowns at the left boundary 
(point L. in FIG 2). 


x L = 0 


( 17 ) 


From equation (9) 

(X L - x A ) = (V A - a A ) (t L - t A ) (18) 

and from equation (10) 

< V L ' Va) ' rir (Pl ‘ ?A> + f V 2 D K aI <*L - *A> = 0 (19) 


With the conditions at A and the total pressure supplied by the reservoir 
known, the pressure (P^), velocity (Vj J, and time (t^) can be found. 

As before, the density was found by integrating the equation of 
state from point A to point L. 

P L = P A ( x + ^7^) (20) 

Boundary Conditions at the Valve 


With reference to FIG 2, the total velocity of the fluid at the valve 
(point. R) is a function of the pressure drop across the valve, and of the 
motion of the vibrating valve. The orifice flow equation, with a 
correction for upstream area, is used to calculate the velocity as a 
function of the pressure drop across the valve (Ref. 18), 

V R = f 2g (p R- 

' A p r (i - 

Because the valve was excited sinusoidally by the vibration exciter, 
the velocity oscillations induced in the water were directly proportional 
to the velocity of the vibration exciter and to the ratio of the projected 
area of the vibrating valve to the area of the pipe line. The displacement 
of the exciter is 


R) ‘ 
5a 

A 2 / 


1/2 


( 21 ) 


X = - ~ cos (2irFt + c|> ) 


( 22 ) 
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The exciter .velocity is found ,as 


— = it dF sin (2^17 Ft + cb) 
dt ‘ 


( 23 ) 


The ratio of the projected area of the valve to the area of the suction 
Therefore, the velocity oscillation caused by 




line 
vibration is 


V 


-R. 


■('* x) 


TrdF sin (2irFt + <j> ) 


(24) 


The total fluid velocity at the valve face is then given as 


Vr = V R] +Vr 2 =-j*. 


Zg (Pr . .R) 

P Rd -0t) J 

A 2 J 


-i 1 /2 


(‘ ■ t) 


(25) 

irdF sin (2 -nr Ft + <j> ) 


From equations (7) and (8) for the C + characteristic 
X R ' X M 


*R - % + 


CVm + a M> 


V R ‘ V M + p^a M (P R " P M> + M 2D ^ (t R ‘ = 0 

also at the valve face 


(26) 

(27) 


-x R = l (28) 

The density relationship is again found by integrating the equation of 
state, this time between and x-^ 

P R = P M (l + PR ' PM ) < Z9 > 

The simultaneous solution of these equations will then yield Vr, 

tj^ and 
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Boundary Conditions for Branched and Compou nd Bines 

Surge pressures in other piping systems can be computed using 
the characteristic equations (7), ( 8 ), (9) and ( 10 ) for each uniform 
section, with the boundary conditions for the case involved. For 
example, the boundary conditions in FIG 3 require that the flow into 
an intersection is the same as the flow leaving the intersection, and 
that the pressure is constant across the intersection. The equations 
of the boundary conditions would then be: 


At junction A V, D* = V 2 D 2 2 (30) 

P iA =P 2 A (31) 

At junction B 

V 2 D 2 2 = V 3 D 3 2 + V 4 D 4 2 (32) 

P 2B = P 3 B = Pj4 B (33) 


Initial Flow Conditions in the Pipeline 


Referring to. FIG 2 before calculations can proceed, a set of 
points must be determined where all of the variables P, V, p, x, and 
t are known. The points used in the analysis are the values when t=0, 
the steady state flow conditions just as valve closure begins. Because 
the conditions were steady state = 0 ), the wave equations ( 2 ) and 
( 6 ) reduce to: 


VfY + f vlvi = 0 

dx 2D p dx 




dP 

Combining the two equations to eliminate yields 

dx 


dV 

-V 


g(3 cLV f 

pV V 2 2D 


dx = 0 


(34) 

(35) 


(36) 
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Because the flow is steady and the suction line area is constant, from 
the requirement of continuity 


G = pV = CONSTANT 

il.iiw + ldx = o 

V G V 2 2D 


Integrating this equation from x to l 


r Y °,l 

/ dV 

Ar 


V 


,v Q i 

SJL/ dV 


Vo, x V 2 2D„ X 


7 dx = 


V Q , x 

The notation V 0f x denotes the velocity at time zero, position x. 


In 


XqJ + M./JL I 

Vo,x G \y a j V o. 


t)' 


2D 


U - x) 


(37) 

(38) 


(39) 


(40) 


Therefore using the above equation and the known conditions at x = / 
the velocity at any point along the .suction line may be calculated. 


'The pressure at any point along the suction line may be found by 
integrating equation (35) 


dP = - 


(41) 

r p o,l 

/ dP = - 



o, X 

“Vq.x V 

(42) 

P 6 - P = 

0,1 o, X 

-p in V °d 
V o.. x 

(43) 


And since' the flowrate is constant, the density variation along the 
pipeline may be found 

G =Po, X V 0)X (44) 

This then defines the flow under steady state conditions. 
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Steps of Computations' 


In solving the characteristic equations, a matrix form was used 
to keep the variables orderly (FIG 4). Here (a) is a variable analagous 
to time, and (b) is a variable analagous to the location along the pipe- 
line, Subscripts are used to denote the location of the variables within 
the matrix. For example P, 2 means the pressure at a = 1, b = 2. 

The steps in computing the unknown quantities are: 

a. Divide the suction line into sections. The number (n) of sec- 
tions chosen will of course determine the degree of accuracy. 

b. Set t = 0 for a--= 0 and solve the steady state equations (40), 

(43) and (44) for b = 0 through n. This sets up the initial conditions. 

c. The unknown quantities P, V, p, and t are then determined 
for the matrix point (1, 0) using the reservoir boundary conditions 
(equations 16 and 17) and the C" characteristic passing through the 
matrix point (0, 1) (equations 18 and 19). Then density at point (1,0) 
is found from equation (20). 

d. The unknowns are then determined for a=l, b~2 through (n-1) 
by using both characteristic curves and the known values at assO. For 
example,, the variables are found at matrix point (1, 3) by simultaneous 
solution of the C + characteristic passing through the matrix point (0, 2) 
arid the C" characteristic (equations 9 and 10) passing through the 
matrix point (0, 4). 

e. The variables are then determined for the matrix point (1, n) 
using the valve boundary conditions (equations 25 and 28) and the C* 
characteristic (equations 7 and 8) passing through the matrix point 
(0, n-1). The density at point (1 3 n) is found from equation (29). 

f. - Using the results calculated for a = 1, the variables for a = 2 
are determined (steps c through e). The calculation then proceeds in 
this manner as far as desired. 

Limitations of Theory 

The limitations of the analysis, as applied to the experimental 
facility, are: 

a. The test pipeline contains a bellows that was not considered in 
the analysis. 
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b. The analysis assumes an abrupt area change between the pipe- 
line and the reservior so that 100% wave reflection is obtained. The 
experimental facility did not have an abrupt area change. 

c. The analysis assumes only the friction losses in the shear 
stress between the fluid and the pipeline, omitting the internal viscous 
dispersion within the fluid. 

d. Two phase fluids are not considered so that once vapor pres- 
sure' is reached, or a bubble exists within the fluid, the equations are 
invalid. 


EXPERIMENTAL INVESTIGATION 
Description of Experimental Facility 

The experimental program was established to measure the tran- 
sient pressure response when a vibrating valve is closed. The 
Experimental Facility is shown in FIG 5. This facility had a uniform 
pipeline of 4- inch O. D. x . 0625 -inch wall aluminum tubing attached 
to the bottom of a 4950 gallon reservoir that could be pressurized to 
20 p-sig. At the lower end of the pipeline, 276 inches below the 
reservoir, a pneumatically operated 4-inch butterfly valve and a 4- 
inch diameter bellows were attached. The butterfly valve was bolted 
to a 22, 000 pound force vibration exciter. A 6-inch drain line with a 
vent was attached to the bottom of the butterfly valve using a bellows 
to allow the butterfly valve to be vibrated as it was closed. 

A second pipeline configuration (see FIG 6) was also used. This 
pipeline was a 6-inch O. D. x 3/l6~inch wall aluminum tube that 
extended from the reservoir halfway to the butterfly valve. The pipe- 
line then was reduced from 6-inch O. D. to 4-inch O. D; the remaining 
length of pipeline to the valve was 4 -inch O. D. tubing. 

Instrumentation 


Four pressure 'transducers were located along the pipeline-to 
collect data. Accelerometers were mounted on the vibration exciter 
and on the butterfly valve. A position indicator was used to measure 
the angular position of the butterfly in the valve. These measurements 
were recorded on an oscillograph using a paper speed of 10 inches/ sec. 
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Also measured during the tests were water temperature, ullage pres- 
sure and liquid level in the reservoir, and the frequency and displace- 
ment of the vibration exciter. 

Control of Test Variables 


The most important factors that determine the maximum surge 
pressure are the time of valve closure, and the variation in fluid 
velocity while the valve is closing. The time of valve closure was 
controlled by orifices on the exhaust port of the pneumatic piston that 
closed the butterfly valve. The variation in fluid velocity while the 
valve was closing could not be controlled., A turbine flowmeter lacks 
sufficient response to measure this transient velocity. Because the 
effective flow area during valve closure could not be directly measured, 
it was approximated by measuring the effective flow area for various 
angles of the butterfly valve during steady flow and assuming the effec- 
tive flow area was the same during unsteady flow. 


To determine the effective flow area during steady flow, a turbine 
flowmeter was installed upstream of the valve. The flowrate through 
the valve and the pressure drop across it were measured for several 
valve positions. The orifice flow (Ref. 18), with the correction for - 
upstream area was used to determine the effective flow area for each 
valve position. 


Q = A t 


1/2 

2g (AP) 1 



solving this for A^. 


A t 


Q 



(45) 


(46) 


The natural frequency of a fluid system is defined as the frequency 
at which a maximum pressure oscillation occurs for a given velocity 
variation of the valve f/ AP \ m ax 1 . The natural frequencies of both 
1 \AVf \ 

pipelines, with the valve closed, were determined experimentally by 
measuring the pressure oscillation at the valve for a given velocity input. 
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The frequencies of vibration of the valve during closure were then 
chosen at and between the natural frequencies. 

With the valve fully- closed, a small velocity oscillation could pro- 
duce a large pressure oscillation, especially at the natural frequencies 
of the systems. Therefore, two levels of vibration were used, The 
first, a low level vibration, was allowed to continue after the valve was 
closed. However, the second, a high level vibration, was stopped when 
the valve fully closed. 

Description of Tests 

Tests were conducted as described below: 

With the butterfly valve closed, the reservoir was filled with water 
and pressurized to. 15 psig. For the low level vibration tests, the 
recorders were started, the vibration exciter was set at the predeter- 
mined frequency and displacement, and the butterfly valve was opened. 
When steady state flow was established, the valve was closed at the 
predetermined rate. The recorders were stopped after the transient 
pressures had damped. The procedure was the same for the high level 
vibration tests except the vibration exciter was turned on after the valve 
was open, and was turned off as soon as the valve was fully closed. 

EXPERIMENTAL RESULTS 
Tests of Uniform 4-Inch Pipeline 


The natural frequencies of the 4-inch pipeline, between 20 and 100 
cps, were determined and are shown in FIG 7. The fundamental natural 
frequency is at 32 cps, with another resonant point at '96 cps. Fre- 
quencies of 32, 48, 64, 80 and 96 cps, were chosen and used as the 
vibration frequency of the valve during closure. 

The steady state effective flow area of the butterfly valve was 
determined and is shown in FIG 8 as a function of the angular position 
of the valve. The effective flow area of the fully open valve is seen to 
be 8. 34 in 2 , compared to a pipeline flow area of 11. 79 in 2 , The steady 
state approximation for a . 250 sec closure was then obtained, as shown 
in FIG 9, by comparing the angular position during closure to the steady 
state effective flow area of the valve, 
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The pressure surges created by valve closures with the 4-inch pipe- 
line are shown in Table I. These maximum pressure surges were 
measured at the pressure transducer located just upstream of the valve. 

A typical pressure response, FIG 10, shows the static pressure at 
the valve for a . 265 sec valve closure with (a) no vibration, (b) low 
level vibration of 32 cps, . 010-inch displacement, and (c) high level 
vibration of 32 cps, . 380 -inch displacement. 

With no vibration, the static pressure increased from the steady 
flow value of 20. 4 psia to a maximum of 117. 1 psia. When the valve 
fully closed, the pressure dropped rapidly and oscillated about the 
steady state level of 48 psia until it was damped out by the losses of 
the system. The pressure surge was determined to be 96.7 psi. At 
the low level vibration of 32 cps, .010 inch dispiaxement, the results 
were similar. The vibrating valve caused only a small- variation in 
the pressure response during the valve closure. However, the continu- 
ation of valve vibration after closure caused a pressure oscillation, 
the first cycle of which was 25. 4 psia to 64. 7 psia. After several 
seconds of vibration* a steady state oscillation as shown in FIG 11 was 
obtained. Although the velocity input to the valve was sinusoidal, the 
pressure oscillation measured at the valve was not sinusoidal. This 
suggests the presence of gas or vapor bubbles in the liquid. 

When the valve was closed during high level vibration of 32 cps, 

. 380 -inch displacement (FIG 10), a significant change in the maximum 
surge pressure was obtained. As the valve closed, the vibration caused 
a pressure oscillation which reached 303. 6 psia, giving a surge of 273. 2 
psi. The minimum pressure during closure was very small (3. 5 psia). 

At the minimum point the curve is flat, probably because of the forma- 
tion of gas bubbles in the fluid. 

The results of a . 250 sec valve closure during vibration at 48' cps 
are shown in FIG 12. The results here are essentially the same as for 
32 cps. The effect of vibrating at a frequency other than the natural 
frequency was twofold; first, it changed the number of pressure oscilla- 
tions during the closure, and secondly, it required vibration at a higher 
level to reach the same magnitude of pressure oscillation. This second 
phenomenon can be seen in FIG 7; as the vibration frequency moves away 
from the natural frequency, a smaller pressure oscillation is obtained 
for the same variation of velocity. At a slower rate of valve closure, 
the pressure response wave lengthened and the maximum pressure surge 
was reduced (Table I). 
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Tests of Pipeline Having an Area Reduction 


The natural frequencies of the pipeline having an area reduction 
were determined and are shown in FIG 13. The fundamental frequency 
is 36. 5 cps, with another resonant point at 82 cps. From this figure, 
frequencies of 30, 36. 5, 48, 59, 70, 82, and 96 cps, were chosen and 
used as the vibration frequency of the valve during closure. 

The pressure surges created by valve closures are shown in 
Table II. A typical pressure response, FIG 14, shows the static pres- 
sure at the valve during a . 232- sec closure with (a) no vibration, (b) 
low level vibration at 70 cps, . 025-inch displacement, and (c) high 
level vibration at 70 cps, . 080-inch displacement. The results here 
are similar to those of the uniform 4- inch suction line, except the 
shape of the pressure response curve was altered slightly. 

ANALYTICAL RESULTS 

The equations describing the transient flow through a pipeline were 
programmed for a digital computer. The boundary conditions of a uni- 
form suction line between a reservoir and a valve were used. The 
equations were then solved using the physical constants of the 4-inch 
pipeline, and the steady state approximation of the effective flow area 
during valve closure (FIG 9). The resulting transient pressure at the 
valve face for a .250 sec closure is shown in FIG 15. Comparing this 
to the experimental results (FIG 10), a higher pressure is predicted 
than measured, and the wave shape is also different. A slight modifi- 
cation of the curve of effective flow area, as shown by the dashed line 
in FIG 9, will yield the correct wave shape and maximum pressure of 
FIG 16. Since this slight modification gave the correct wave shape,, 
the assumption previously made, that the steady state effective flow 
area was equal to the transient effective flow area, is only approxi- 
mate. For comparison, the transient pressure created by an instan- 
taneous closure is predicted [Equation (1)] as 1760 p si. 

FIG 17 shows the calculated results under conditions similar to 
those of FIG-16 except that the valve was vibrating at 32 cps, with a 
displacement of .030 inches. The vibration only slightly modifies the 
base curve. However, when the displacement is increased to . 380 
inches (FIG 18), there is a marked increase in the maximum pressure 
during the valve closure. These results compare favorably with the 
experimental results shown in FIG 10. FIG 18 shows that the pressure 
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decreases quite rapidly, and according to the equations, the pressure 
would continue decreasing to below zero psia. Of course, this could 
not occur, so computation was stopped. This is a significant point, 
because if the vapor pressure of the fluid is reached, the fluid will- 
cavitate, similar to the cavitation in ultrasonic cleaning. This is 
probably what happened in the experimental results shown in FIp 10. 
When cavitation does occur, the expressions in- this analysis do not 
apply because the assumption of a constant bulk modulus for the fluid 
is invalid. 

The pressure variation at the valve for a sinusoidal vibration of 
32 cps with . 030 inches displacement with the valve closed is shown 
in FIG 19. The pressure variation steadily increases, damped only by 
the losses of the system. Without losses or cavitation, the oscillation 
would increase undamped for the natural frequency. For other than 
the natural frequencies, there would be a steady state level even with- 
out losses. If the pressure oscillations caused cavitation, large local 
pressures could be created due to collapse of the bubbles formed. 

CONCLUSIONS 

Analytical calculations, which are verified by results of the 
experimental tests, show that a vibrating valve can cause a substantial 
increase in the maximum transient pressure in a pipeline during valve 
closure. ~The magnitude of the transient pressure is directly related 
to the magnitude of valve vibration, and is also a function of the fre- 
quency of forced vibration. If the valve is closed, small valve vibrations 
at the natural frequency of the fluid system cause large pressure oscil- 
lations. 

The initial fluid velocity, time of valve closure, and the rate of 
decrease in flow area are important factors in determining the maxi- 
mum transient pressure. The magnitude of the transient pressure of 
a fluid system, with or without vibration, may be calculated by the 
methods reported if the boundary conditions of the fluid system are 
known. This method does not neglect the nonlinear terms of the wave 
equations (including friction), and it provides results along the suction 
line as well as at each end. 
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TABLE I 


FIG 10 


FIG 12 


Measured Pressure Surges 
Uniform 4- Inch Pipeline 


Test 

No. 

Valve 

Closure Time 
(sec) 

Vibration 

Frequency 

(cps) 

Vibration 
Di splacement 
(in) 

Pressure Surge 
at Valve Face 
(psi) 

2 

.268 

0 

0 

97' 

14 

.270 

32 

. 010 

97 

74 

. 260 

32 

. 380 

304 

17 

.250 

48 

. 038 

110 

71 

. 250 

48 

. 170 

165 

30 

. 260 

64 

. 042 

101 

76 

. 251 

64 

. 095 

106 

33 

.251 

80 

. 024 

102 

67 

. 265 

80 

. 060 

116 

46 

. 264 

96 

. 010 

102 

66 

. 243 

96 

. 043 

154 

80 

. 504 

0 

0 

52 

102 

.449 

32 

. 010 

59 

120 

.466 

32 

. 380 

160 

103 

. 457 

48 

. 038 

69 

117 

.462 

48 

. 170 

116 

90 

.481 

64 

. 042 

60 
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TABLE I (Concluded) 


Test 

No. 

Valve 

Closure Time 
(sec) 

Vibration 

Frequency- 

Ceps) 

Vibration 

Displacement 

(in) 

Pressure Surge 
at Valve Face 
(psi) 

112 

. 475 

64 

. 095 

70 

92 

. 477 

80 

. 024 

64 

111 

. 473 

80 

. 060 

122 

94 

. 483 

96 

. Oio 

64 

106 

. 460 

96 

. 043 

no 

4 

• 998 

0 

0 

40 

10 

1. 025 

32 

. 010 

56 

49 

. 939 

32 

. 380 

148 

21 

. 955 

48 

. 038 

62 

52 

. 944 

48 

. 170 

75 

24 

. 957 

64 

. 042 

44 

55 

. 948 

64 

. 095 

156 

37 

. 967 

80 

. 024 

51 

59 

. 958 

80 

. 060 

165 

43 

. 948 

96 

. 010 

54 

60 

. 932 

96 

. 043 

113 
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TABLE -II 


Measured Pressure Surge 
Pipeline with Area Transition 


Test No. 

Valve 

Closure Time 
(sec) 

Vibration 

Frequency 

(cps) 

Vibration 

Displacement 

(in) 

Pres sure 
Surge 

(psi) 

233 

. 127 

0 

0 

162 

235 

. 129 

30 

. 030 

166 

305 

. 118 

30 

. 300 

192 

238 

. 130 

36. 5 

.005 

162 

307 

. 118 

36. 5 

.295 

193 

242 

. 130 

48 

. 030 

159 

310 

. 132 

48 

. 150 

181 

244 

. 129 

59 

.045 

1 55 

313 

. 135 

59 

. 100 

165 

248 

. 129 

70 

. 025 

160 

317 

. 131 

70 

. 080 

171 

251 

. 130 

82 

. 004 

155 

319 

. 125 

82 

. 058 

176 

254 

. 129 

96 

. 015 

158 

322 

. 124 

96 

. 042 

175 

230 

. 238 

0 

0 

93 

210 

. 224 

30 

. 030 

95 
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TABL.E II (Continued) 


Test No. 

Valve 

Closure Time 
(sec) 

Vibration 

Frequency 

(cps) 

Vibration 

Displacement 

(in) 

Pressure 

Surge 

(psi) 

325 

. 238 

30 

.280 

120 

212 

. 222 

36. 5 

. 006 

86 

328A 

. 241 

36. 5 

. 250 

139 

215 

.224 

48 

. 030 

90 

332 

. 224 

48 

. 150 

106 

218 

. 223 

59 

. 045 

90 

334 

.239 

59 

. 100 

103 

222A 

.235 

70 

,025 

102 

337 

. 222 

70 

. 080 

126 

223 

. 230 

82 

.004 

93 

341 

.234 

82 

. 060 

150 

228 

. 228* 

96 

. 015 

101 

345 

. 228 

96 

. 042 

149 

278 

. 549 

0 

0 

49 

274 

. 546 

30 

. 030 

55 

348 

. 561 

30 

. 300 

137 

271 

. 550 

36. 5 

. 005 

51 

351 

. 567 

36. 5 

. 250 

190 
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TABLE II (Concluded) 


Test No. 

Valve 

Closure Time 
(sec) 

Vibration 

Frequency 

(cps) 

Vibration 

Displacement 

(in) 

Pressure 

Surge 

(psi) 

270 

. 535 

48 

.030 

68 

352 

. 569 

48 

. 150 

108 

267 

. 548 

59 

. 045 

54 

357 

. 532 

59 

. 100 

127 

264 

. 541 

70 

. 025 

60 

358 

. 559 

70 

.080 

69 

260 

. 543 

82 

.004 

50 

362 

.563 

82 

.060 

81 

257 

.*541 

96 

.015 

78 

364 

. 565 

96 

. 042 

73 

385 

.843 

30 

. 300 

67 

382 

.853 

36. 5 

.250 

72 

380 

. 839 

48 

. 150 

62 

376 

. 821 

59 

. 100 

61 

375 

. 820 

70 

. 080 

70 

371 

. 826 

82 

. 060 

66 
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1. Pressurization Supply Valve 

2. Pressurization Regulator 

3. Vent & Safety Relief Valve 

4. Drain Valve 

5 . Surge Valve 

6. Fill Valve 

7 . Reservoir 

8. Vibration Exciter 

PI - Pressure Below Valve 
P2 - Pipeline Transient Pressur 
P3 - Pipeline Transient Pressur 
P4 - Pipeline Transient Pressur 
P5 - Liquid Level Indicator 
P6 - Reservoir Pressure 
T1 - Fluid Temperature 


Position - 
Indicator 


FIGURE 5. SCHEMATIC OF TEST FACILITY 
(4-INCH DIAMETER PIPELINE) 


R-P&VE-PE 

12-1-64 
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Va lve 
Angular 
Position . 
(% closed) 



DO 

NO 


FIGURE 8. STEADY STATE EFFECTIVE FLOW ARE'A OF BUTTERFLY VALVE 
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FIGURE 12. MEASURED TRANS IEKT PRESSURE, , 250 SEC CLOSURE, ‘48 CPS 
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TRANSIENT PRESSURE, PIPELINE WITH AREA TRANSITION 
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FIGURE 18. CALCUIATED PRESSURE, ,250 SEC CLOSURE, 32 CPS, .380 INCH 





APPENDIX A 


DERIVATION OF THE MOMENTUM EQUATION 
[Equation (2)] 

Consider fhe .control -volume of length dx in the uniform pipeline 
of FIG 1. At. face a of the control volume, the pressure (P) acts 
over the area (A,) exerting a force (PA) on .the control volume. 


F a = PA 


(47) 


>‘(The direction of positive force is the same as the direction 
•of .positive velocity) 

gjp 

A.t face b, the pressure (P + dx) acts over the area (A) exert- 

ing a force on -the -control volume of 

3P 

F b =_(p + sr dx > J < 48 > 

A shear force is also exerted on the control volume by .the 
■fluid moving pa-st the pipe wails. The direction of -the shear force 
i<s opposite «the direction of -.the fluid velocity, and its magnitude is 
A w t w . {(W'he-re A w is .the- wetted area of the control volume and Toj 
is the' .shearing stress exerted on the stream by the wails). The 
wetted *area -expressed in terms of the pipeline dimensions is 

_ , /4D\ /tt D\ 4dx 4A 

A u = nDdx^ j = [—) -Q- = - 5 " dx (49) 

The Darcy- Weisbach friction factor is defined as four times the rati ; o 
of the '.wall .’shearing stress to the dynamic head of the fluid 


4 

2 /2g 


(50) 


The magnitude of the shear force in terms of the friction factor 


4A dx f pY 2 
D 2g4 


(51.) 
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The direction of the shear force is controlled by using the fluid 
velocity times its absolute magnitude (VIVI) instead of V 2 . (VIVI). 
always has the sign of V, whereas V 2 is always positive. The shear 
force acting on the control volume then becomes 


F = T dx fpVjvl 

s D 2g 


(52) 


Summing the forces on the control volume 

ZF = F a + F b + F s = AP - A (P + -|£ dx) - A ^ fp ^ V l 

The mass of the control volume (which is constant) is 

M = _JL Adx 
g 

The acceleration of the control volume is 



(53) 


(54) 


(55) 


In unsteady flow, the velocity is a function of both the location along 
the pipeline (x) and the time (t) 




V = 

f(x, t) 

(56) 

dV = 

av 

dx + . 

31 dt 

(57) 


dx 


at 

dV _ 

av 

dx 

av 

(58) 

dt 

3x 

>dT * 

at 

dx = 

Y 



(59) 

dt 




_ dV 


av 

av 

(60) 

dt 

= Y ■ 

-5 — + 
9x 

at 
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Using Newton’ s second 5 law 


tAP* - A (P- + 


pvlvl 

D 2g 


Simplifying 


„ 9V 9V g 3P f , , 

V + ~ + f S? + 2D V|V| 
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APPENDIX B 


DERIVATION OF THE CHARACTERISTIC EQUATIONS 
[Equations (7), (8), (9), & (10)1 


The basic equations used are the wave equations: 

M + y *L + jl _|P + 4 vlvi =0 

9t 9x p 9.x 2D 


9P 


+ P 

9t 1 9 x 


3V + V |P = o 


3x 


(2) 

( 6 ) 


A linear combination of these equations yields 


+V JY + 1 ± V |v|+«— + *R - 9V 

9x P 9x 2D 9t 

Rearranging this equation yields 


- + aV = 0 

ox 9x 


+ (v + ) 

9V 1 . 

+ O' 

+ (V + 

-JUihT 

j_8t 

9x J 

l 9t 

cap 9x J 


2D 


(63) 


(64) 


Since P and V are the dependent variables, and x and t are the 
independent variables. 


dP 9P dx t 9P 

dt 9t + ‘dt '9x 

and 

dV 9V dx 9V 

dt 9t dt 9x 


(65) 


( 66 ) 


In equation (64), the total derivations and can replace the 

quantities in brackets if 


±5=V+G'p=V + -L (67) 

dt cap 
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This requires that 


* • i(*-) 


1/2 


( 68 ) 


Therefore, two equations were obtained, one for the positives 
(defining the C + characteristic) and .another for the negative or 
(defining the C~ characteristic). 


Substituting the positive a into equation (67) yields 

, / \ .1/2 
~ = y+ /-M.) 

Since the velocity of wave propagation is defined as 

/ \ 1 /2 


(69) 


(70) 


.equation (69) becomes 


dx 

'diT 


V’+ a 


(71) 


Substituting this into the total differential equations (65) and (66) 


dP -3P , 9P 

— — = •— + (V + a) — - 

dt at 8x 


= JV + (V + a) 4 2L 

dt 8t 3x 


(72;) 

(73.) 


Equation (64) with the positive a substitution becomes 


■ ~ H 
|_ a - 


™ + ( V + a) iX 
ft ax 


. + — i- _ 

pa _ 


t (V + a) ^ 

at ax 


_ + 2 3 v|y| = 0 


(74) 
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substituting (72) and (73) into equation (74) 


dV + dP + £ 

dt pa dt 2D 


V |v I = 0 


therefore along C + 


dx = (V + a) dt 


dV + 


pa 


dP + 25 


V|v|dt = 0 


similarly along C" 

dx = (Y - a) dt 

d v - -fr dP + vivrdt = 0 


(75) 

(76) 

(77) 

(78) 

(79) 


These characteristics are shown in FIG 2. If the velocity and pressure 
are known at any two points on the x-t diagram (points A and B for 
example) the characteristics can be solved for the conditions at point 
o. To do this, a linear approximation- between A and o and B and o, 
as shown below, is used. 



dx = f (x^ (x 2 - x 1 ) 


The characteristics then become 


.(80) 


Along C* 

( x o ~ X A> = < V A + a A } (t o “ l A> ( ? ) 

g f i j 

(V o - V A ) + 7^2 {P o “ P A) + 2D Y a |v a | (t Q - t A ) = o 

( 8 ) 
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Along. C~ 


(x o - x b) = < V B - a B> (‘o - %) (9) 

< V ° - V B ) - (Po - P B > + 2 3 V bKb| (to ' t B )=0 (10) 

With the conditions given at points A and B, these equations may be 
solved simultaneously for P 0 , V 0> x Q and t Q . 
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